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Abstract 

Motivated by the increasing evidence for the need of a geometry that resembles Bianchi morphology 
to explain the observed anisotropy in the WMAP data, we have discussed some features of the Bianchi 
type-y/o universes in the presence of a fluid that wields an anisotropic equation of state (EoS) parameter 
in general relativity. We present two accelerating dark energy (DE) models with an anisotropic fluid in 
Bianchi type-V^/o space-time. To prevail the deterministic solution we choose the scale factor a(t) = \/t"e', 
which yields a time-dependent deceleration parameter (DP), representing a class of models which generate 
a transition of the universe from the early decelerating phase to the recent accelerating phase. Under the 
suitable condition, the anisotropic models approach to isotropic scenario. The EoS for dark energy u is 
found to be time-dependent and its existing range for derived models is in good agreement with the recent 
observations of SNe la data (Knop et al. in Astrophys. J. 598:102, 2003), SNe la data with CMBR anisotropy 
and galaxy clustering statistics (Tegmark et al. in Astrophys. J. 606:702, 2004) and latest combination of 
cosmological datasets coming from CMB anisotropics, luminosity distances of high redshift type la supernovae 
and galaxy clustering (Hinshaw et al. in Astrophys. J. Suppl. 180:225, 2009, Komatsu et al. in Astrophys. J. 
Suppl. 180:330, 2009). For different values of n, we can generate a class of physically viable DE models. The 
cosmological constant A is found to be a positive decreasing function of time and it approaches to a small 
positive value at late time (i.e. the present epoch) which is corroborated by results from recent type la 
supernovae observations. We also observe that our solutions are stable. The physical and geometric aspects 
of both the models are also discussed in detail. 
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1 INTRODUCTION 

Recent cosmological observations obtained by SNe la (Riess et al. 1998; Perlmutter et al. 1999) suggested that 
the expansion of the universe is accelerating. Recent observations of SNe la of high confidence level (Tonry et 
al. 2003; Riess et al. 2004; Clocchiatti et al. 2006) have further confirmed this. In addition, measurements of 
the cosmic microwave background (CMB) anisotropics (Bennett et al. 2003; de Bernardis et al. 2000, Hanany 
et al. 2000), large scale structure (LSS) (Tegmark et al. 2004a, b; Spergel et al. 2003), the Sloan digital sky 
survey (SDSS) (Seljak et al 2005; Adelman-McCarthy et al. 2006), the Wilkinson microwave anisotropy probe 
(WMAP) (Perlmutter et al. 2003) and the Chandra x-ray observatory (Allen et al. 2004) strongly indicate that 
our universe is dominated by a component with negative pressure, dubbed as dark energy, which constitutes 
with ~ 3/4 of the critical density. The cosmic acceleration is realized with negative pressure and positive energy 
density that violate the strong energy condition. This violation gives a reverse gravitational effect. Due to this 
effect, the Universe gets a jerk and the transition from the earlier deceleration phase to the recent acceleration 
phase take place (Caldwell et al. 2006). A recent survey of more than 200, 000 galaxies appears to confirm the 
existence of dark energy, although the exact physics behind it remains unknown ( Rincon 2011). 

During the last two decades cosmology is speedily becoming an experimental involvement of physics. The 
theoretical models can be tested, and new and more accurate data in the near future will restrict our conceptions 
of the Universe to within few percent accuracy. The simplest candidate for the dark energy is the cosmological 
constant (Overduin and Cooperstock 1998; Sahni and Starobinsky 2000; Komatsu et al. 2009; Kachru et al. 
2003) which suffers from conceptual problems such as fine-tuning and coincidence problems (Weinberg 1989). 



Other scenarios include, Quintessence (Wetterich 1988; Ratra and Peebles 1988), Chameleon (Khoury and Welt- 
man 2004), K-essence (Chiba et al. 2000; Armendariz-Picon et al. 2000), which is based on earlier work of 
K-inflation (Armendariz-Picon et al. 1999), modified gravity (Capozziello 2002; Caroll et al. 2004; Nojiri and 
Odintsov 2003, 2004; Abdalla et al. 2005; El-Nabulsi 20011a), Tachyon (Padmanabhan 2002) arising in string 
theory (Sen 2002), Quintessential inflation (Peebles and Vilenkin 1999), Chaplygin gas as well as generalized 
Chaplygin gas (Srivastava 2005; Bertolami et al. 2004; Bento et al. 2002; BiUc et al. 2002; AveUno et al. 2003), 
cosmological nuclear energy (Gupta & Pradhan 2010). Recently, El-Nabulsi (2011b), Feng and Yang (2011), 
Biesiada et al. (2011), Singh and Chaubey (2012), Amirhashchi et al. (2011a) and Pradhan & Amirhashchi 
(2011a) have studied DE models in different context. In spite of these attempts, still cosmic acceleration is a 
challenge to modern cosmology and modern astrophysics. 

In general relativity, the evolution of the expansion rate is parameterized by the cosmological equation of 

state (the relationship between temperature, pressure, and combined matter, energy, and vacuum energy den- 
sity for any region of space). Measuring the equation of state for dark energy is one of the biggest efforts in 
observational cosmology today. The DE model has been characterized in a conventional manner by the equation 
of state (EoS) parameter Lo{t) = ^ which is not necessarily constant, where p is the energy density and p is the 
fluid pressure (Carroll and Hoffman 2003). The present data seem to slightly favour an evolving dark energy 
with EoS a; < — 1 around the present epoch and w > — 1 in the near past. Obviously, w cannot cross —1 for 
quintessence or phantom alone. Some efforts have been made to build a dark energy model whose EoS can 
cross the phantom divide. The simplest DE candidate is the vacuum energy {ui = —1), which is mathemati- 
cally equivalent to the cosmological constant (A). The other conventional alternatives, which can be described 
by minimally coupled scalar fields, are quintessence (cj > — 1) (Steinhardt and Wesley 2009), phantom energy 
(w < — 1 (Caldwell 2002) and quintom (that can across from phantom region to quintessence region as evolved) 
and have time dependent EoS parameter. Some other limits obtained from observational results coming from 
SNe la data (Knop et al. 2003) and combination of SNe la data with CMBR anisotropy and galaxy cluster- 
ing statistics (Tegmark et al. 2004a,b) are —1.67 < ui < —0.62 and —1.33 < ui < —0.79, respectively. The 
latest results in 2009, obtained after a combination of cosmological datasets coming from CMB anisotropies, 
luminosity distances of high redshift type la supernovae and galaxy clustering, constrain the dark energy EoS 
to -1.44 < LO < -0.92 at 68% confidence level (Hinshaw et al. 2009; Komatsu et al. 2009). However, it 
is not at all obligatory to use a constant value of oj. Due to lack of the observational evidence in making a 
distinction between constant and variable w, usually the equation of state parameter is considered as a constant 
(Kujat et al. 2002; Bartelmann et al. 2005; Yadav 2011) with phase wise value — 1,0, — | and -|-1 for vacuum 
fluid, dust fluid, radiation and stiff dominated universe, respectively. But in general, w is a function of time or 
redshift z or scale factor a as well (Ratra and Peebles 1988; Jimenez 2003; Das et al. 2005). In earlier various 
form of time dependent w have been used for variable A models by Mukhopadhyay et al. (2008). Recently, 
dark energy models with variable EoS parameter have been studied by Ray et al. (2010), Akarsu and Kilinc 
(2010a,b), Yadav et al. (2010), Yadav and Yadav (2010), Pradhan and Amirhashchi (2011a), Pradhan et al. 
(2011), Amirhashchi et al. (2011a,b) and Saha & Yadav (2012). In well-known reviews on modified gravity (No- 
jiri and Odintsov 2007, 2011), it is clearly indicated that any modified gravity may be represented as effective 
fluid with time dependent w. The dark energy universe EoS with inhomogencous, Hubble parameter dependent 
term is considered by Nojiri and Odintsov (2005). Later on, Nojiri and Odintsov (2006) have also presented 
the late-time cosmological consequences of dark energy with time-dependent periodic EoS in oscillating universe. 

Today there is considerable evidence, which suggests that the universe may be isotropic and homogeneous. 
After discovery of CMB radiation, cosmology became a precision science. The CMB radiation is also considered 
to be a major experimental evidence on which the most commonly accepted theory about the origin of universe, 
i.e. "Big-Bang" cosmology, is based. Statistical Isotropy (SI) is usually assumed in almost all CMB studies. But, 
now, there exist many indications which suggest that CMB may violate this assumption. Apart from the CMB 
there are some other indications of violation of SI which suggest the existence of a preferred direction in the 
universe. These indications include distributions of polarizations from radio galaxies (Birch 1982; Jain & Ral- 
ston 1999; Jain et al. 2004) and statistics of optical polarizations from quasars (Hutsemekers 1998; Hutsemekers 
& Lamy 2001; Jain et al. 2004; Ralston & Jain 2004). Polarization of electromagnetic waves coming from 
distant Radio Galaxies and Quasars measured at radio and optical frequencies respectively are not consistent 
with the assumptions of SI, rather radio polarizations are organized coherently over the dome of the sky and 
optical polarizations are aligned in a preferential direction on very large scales, violating the assumed isotropy 
of the universe. These study confirmed strong significance of anisotropy and also claimed that the statistics are 
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not consistent with isotropy at 99.9% confidence level. It has also been observed that the quadrupole and the 
octopole have almost all their power perpendicular to a common axis in space pointing towards Virgo cluster 
(Tegmark et al. 2003; de Oliveira-Costa 2004). The dipole, which is commonly attributed to our motion relative 
to the CMB rest frame, also aligns in the same direction as quadrupole and octopole which is not expected under 
the condition of statistical isotropy. Another indication of anisotropy in CMB data is the presence of a cold spot 
with improbably low temperature. It was found by Cruz et al. (2005) by using Spherical Mexican Hat Wavelet 
analysis on WMAP data. Several authors have also searched for anisotropy using the supernova data set. Jain 
et al. (2007) found violation of isotropy in this data. Subsequently, there have been a large number of studies 
(Bielewicz et al. 2004; Eriksen et al. 2004; Katz & Weeks 2004; Bielewicz et al. 2005; Prunet et al. 2005; Bernui 
et al. 2006; de Oliveira-Costa & Tegmark 2006; Freeman et al. 2006; Bernui et al. 2007; Land & Magucijo 2007) 
which claim the CMB is not consistent with isotropy. The possible violation of SI in CMB has lead to many 
theoretical studies. Several physical explanations for the observed anisotropy have been put forward (Cline et al. 
2003; Contaldi et al. 2003; Kesden et al. 2003; Armendariz-Picon 2004; Berera et al. 2004; Gordon et al. 2005; 
Abramo et al. 2006; Campanelli et al. 2007; Rodrigues 2008). Land & Magueijo (2005) found evidence that the 
detected anisotropy has positive mirror parity. The generation and evolution of primordial perturbations in an 
anisotropic universe has also been studied (Armendariz-Picon 2006; Battye & Moss 2006; Pereira et al. 2007) 
along with the possibility of anisotropic inflation (Hunt & Sarkar 2004; Buniy et al. 2006; Donoghue et al. 2009). 

The possible violation of global isotropy in the CMB has been a subject of intense research after the publica- 
tion of WMAP data. The possible alignment of axes corresponding to several diverse data sets in the direction 
of the Virgo cluster makes this extremely interesting. In recent years, there have been a large number of studies 
which claim the CMB temperature fluctuations are not consistent with statistical isotropy and thus questioning 
the cosmological principle. The CMB is considered to be a major experimental evidence which supports the cur- 
rent/present models of the observed universe and from this CMB observations several people found significant of 
anisotropic scenario. Based on these studies one may not preclude the possibility that our universe is anisotropic. 

There is a general agreement among cosmologists that cosmic microwave's background anisotropy in the small 
angle scale holds the key to the formation of the discrete structure. The theoretical argument (Misner 1968) and 
the modern experimental data support the existence of an anisotropic phase, which turns into an isotropic one. 
The anisotropy of the DE within the framework of Bianchi type space-times is found to be useful in generating 
arbitrary ellipsoidality of the Universe, and to fine tune the observed CMBR anisotropics. Koivisto and Mota 
(2008a, 2008b) have investigated cosmological models with anisotropic EoS and have also shown that the present 
SN la data allows large anisotropy. Recently, Akarsu and Kilinc (2010c) have described some features of the 
Bianchi type-I universes in presence of fluid that wields an anisotropic EoS. Hence, for a realistic cosmological 
model one should consider spatially homogeneous and anisotropic space-times and then show whether they can 
evolve to the observed amount of homogeneity and isotropy. The only spatially homogeneous but anisotropic 
models other than Bianchi type models are the Kantowski-Sachs locally symmetric family. See Ellis & van Elst 
(1999) for generalized, particularly anisotropic, cosmological models and Ellis (2006) for a concise review on 
Bianchi type models. The motivation for this investigation comes from the hints of statistical anisotropy of our 
Universe that several observations seem to suggest. 

Bianchi type-V/o (B-VIq) space-time in connection with massive strings is studied by Pradhan and Bali 
(2008) and Bali et al. (2008). Belinchon (2009) studied several cosmological models with B-V/q & HI symmetries 
under the self similar approach. Given the growing interest of cosmologists, here, we propose to study the 
evolution of the Universe within the framework of a B-VIq space-time. Recently, Amirhashchi et al. (2011c) and 
Pradhan et al. (2012) presented dark energy models in an anisotropic B-VIq space-time by considering constant 
and variable deceleration parameter respectively. In this paper, we have investigated two new B-VIq DE models 
with variable co by assuming different scale factors in such a way that they provide time dependent deceleration 
parameter in presence of anisotropic fluid. The out line of the paper is as follows: In Sect. 2, the metric and 
the field equations arc described. Section 3 deals with the solutions of the field equations. Section 4 deals with 
physical and geometric behaviour of the model. Section 5 deals with the stability of the corresponding solutions. 
In Section 6, wc describe an other dark energy model and its physical aspects. In Sect. 7, we again examine the 
stability of corresponding solutions for second DE model. Finally, conclusions are summarized in the last Sect. 
8. 
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2 THE METRIC AND FIELD EQUATIONS 

We consider totally anisotropic Bianchi type-V^/g line element, given by 

ds^ = -df + AHx" + B^e^-dy^ + C^e^^^dz^ (1) 

where the metric potentials A, B and C are functions of t alone. This ensures that the model is spatially 
homogeneous. 

The simplest generalization of EoS parameter of perfect fluid may be to determine the EoS parameter 
separately on each spatial axis by preserving the diagonal form of the energy momentum tensor in a consistence 
way with the considered metric. Therefore, the energy momentum tensor of fluid can be written, most generally, 
in anisotropic diagonal form as follows: 

=dtag[TS,T^,Ti,Tl] (2) 

Allowing for anisotropy in the pressure of the fluid, and thus in its EoS parameter, gives rise to new possibilities 
for the evolution of the energy source. To see this, we first parametrize the energy momentum tensor given in 
([2|) as follows: 

T/ = diag[p, -p.^, -py, -p^] 

= diag[l, -u^, -ujy, -LOz]p 

= dmg[l,-uj,~{iv + 6),-iuj + -f)]p. (3) 

Here p is the proper energy density, Px,Py and pz are the pressures, and ujx^^y and uJz are the directional EoS 
parameters along the x, y and z axes, respectively; uj is the deviation-free EoS parameter of the fluid. The 
deviation from isotropy is parametrized by setting uJx — ^ and then introducing skewness parameters 6 and 7 
which are the deviations from to, respectively along the y and z. u, 5 and 7 are not necessarily constants and 
might be function of the cosmic time, t. 

The Einstein's field equations (with gravitational units, SvrG = 1 and c = 1) read as 

Ri~\R9l^-Ti, (4) 

where the symbols have their usual meaning. In a comoving co-ordinate system, Einstein's field equation (|^. 
with ^ for B-t//o metric ([T]) subsequently lead to the following system of equations: 

~ujp, (5) 

> + 5)p, (6) 

' + 1)P, (7) 

^_ - P, (8) 

Here and in what follows an over dot denotes ordinary differentiation with respect to t. 
The spatial volume for the model ([1]) is given by 

= ABC. (10) 

We define a = (ABC) 3 as the average scale factor so that the Bubble's parameter is anisotropic and may be 
defined as 
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The deceleration parameter q, the scalar expansion 0, shear scalar cr^ and the average anisotropy parameter A, 
are defined by 




where AHi = Hi — H(i = x, y, z). 



3 SOLUTIONS OF THE FIELD EQUATIONS 

Integrating Eq. we obtain 

C^iB, (16) 

where £ is an integrating constant. Now if we put the value of Eq. ([T5|) in ([7|) and subtract the result from Eq. 
(pl), we obtain that the skewness parameters along y and z axes are equal, i.e 6 — "f. 

Therefore, equations dS])-® are reduced to 

A B Ab 1 , , 
I + B + l5-A^ = -^" + ^)^' 
Ab B^ 1 

'ab^b^-^-^- 

The field equations ([T7 | - p^ are a system of three linearly independent equations with five unknown parameters 
A, B, uj, p and 7. Two additional constraints relating these parameters are required to obtain explicit solutions 
of the system. 

In literature it is common to use a constant deceleration parameter (Akarsu and Kilinc 2010a, 2010b; 
Amirhashchi et al. 2011c; Pradhan et al. 2011; Kumar and Yadav 2011; Yadav 2011), as it duly gives a 
power law for metric function or corresponding quantity. The motivation to choose such time dependent DP is 
behind the fact that the universe is accelerated expansion at present as observed in recent observations of Type 
la supernova (Riess et al. 1998; Perlmutter et al. 1999; Tonry et al. 2003; Riess et al. 2004; Clocchiatti et al. 
2006) and CMB anisotropics (Bennett et al. 2003; de Bernardis et al. 2000; Hanany et al. 2000) and decelerated 
expansion in the past. Also, the transition redshift from deceleration expansion to accelerated expansion is 
about 0.5. Now for a Universe which was decelerating in past and accelerating at the present time, the DP must 
show signature flipping (see the Refs. Padmanabhan and Roychowdhury (2003), Amendola (2003), Riess et al. 
(2001)). So, in general, the DP is not a constant but time variable. This motivates to choose such a scale factor 
which yields a time-dependent DP. At this juncture it should be stated that some authors first choose the scale 
factors in power law, exponential or in other form and then find out other variables with some conditions under 
these solutions. 



In this paper, following Saha et al. (2011) and Pradhan & Amirhashchi (2011b), we take following ansatz 
for the scale factor, where increase in term of time evolution is 

a{t) = VFe^, (20) 

where n is a positive constant. Saha et al. (2011) and Pradhan & Amirhashchi (2011b) examined the relation 
(|20p in studying two-fluid scenario for dark energy models in an FRW universe and accelerating DE models 
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in Bianchi type-V space-times respectively. This ansatz generalized the one proposed by Amirhashchi et al. 
(2011b). If we put n = in Eq. (PO)) . it is reduced to a{t) = -s/e* i.e. a exponential law of variation of scale 
factor. This choice of scale factor yields a time-dependent deceleration parameter (see Eq. (pO|) such that 
before DE era, the corresponding solution gives inflation and radiation/matter dominance era with subsequent 
transition from deceleration to acceleration. Thus, our choice of scale factor is physically acceptable. 

It is worth mention here that one can also select many to many ansatz other than Eq. (j20p which mimic 
accelerating universe but one should also be careful to check the physical acceptability and stability of their 
corresponding solutions otherwise does not prove any relation of such solutions with observable universe. Eq. 
([20|) yields physically plausible solutions. 



Secondly, we assume that the expansion {6) is proportional to shear (ct). This condition and Eq. (jl6p lead 

to 

which yields to 

4 - (22) 

where m = ^"^"^fva^ 

ao are arbitrary constants. Above equation, after integration, reduces to 

A = /3(B)", (23) 

where /3 is an integrating constant. Here, for simplicity and without any loss of generality, we assume /3 = 1. 
Hence we have 

A = iB)"\ (24) 

Collins et al. (1980) have pointed out that for spatially homogeneous metric, the normal congruence to the 
homogeneous expansion satisfies that the condition ^ is constant. 

Using equations ([T6| . (f20|) and ([24|) in (fTT|) . we obtain the expressions for metric functions as follows 

B{t) =£i{t"e')'^(^, (25) 
C(t) =^2(i"e*)5(^, (26) 
A{t) ^ esife')^'^ , (27) 

1 

where, £i ~ k , £2 ~ Mi, £3 — £™ and k is an integrating constant. 

Hence the model H]) reduces to 

ds^ ^-dt^+ £l e* ) tSvi) dx^ +£l{re') dy^ 

+ £l{f'e')T^dz^. (28) 

4 PHYSICAL ASPECTS OF THE DARK ENERGY MODEL 

The expressions for the Hubble parameter (H), scalar of expansion (9), shear scalar (cr), the spatial volume {V) 
and the average anisotropy parameter {Am) for the model (j28p are given by 

e = 3H^l(l+'4) , (29) 

(30) 



2 V t 
2n 



^ (n + i)2 

-^4f^V(l + T) > (31) 



4 V m + 2 / V t 
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Figure 1: The plot of energy density p versus t. Here Iq = 0.1, m — 1 



V = {ee')K (32) 
A„,^2(^^)\ (33) 



From Eqs. (P ^ — it is observed that at i = 0, the spatial volume vanishes and other parameters 6, a, 
H diverge. Hence the model starts with a big bang singularity at t = 0. This is a Point Type singular- 
ity (MacCallum 1971) since directional scale factor A{t), B{t) and C{t) vanish at initial time. Since ^ ^ 
except m = 1, hence the model is anisotropic for all values of m except for m ^ \. The dynamics of the 
mean anisotropy parameter depends on the value of to. We observe that when m = 1, = (i.e. the case of 
isotropy). Thus, the observed isotropy of the model can be achieved in cosmological constant region (see, Fig. 2). 

The energy density of the fluid can be find by using Eqs. ((25|) & p7)) in 



where = 77- Using Eqs. (pSj). (P71) and ([M)) in pT)) . the EoS parameter w is obtained as 



jjirw (1 + t)' - (;ot + e^)'^ 
— TTT — rrr~ 79 ■ y"^^) 



£o(t«e*)-T^ + |(|r^)( 



1 + ^ 



Using Eqs. (1^ . (P7|) . and ([55)1 in p^. the skewness parameters 5 (or 7) (i.e. deviations from w along 
y-axis and z-axis) are computed as 

!(^){(i + f)'-i?}-24(tv)-™ 

= 7 = ^ ^-^^ ^ /„ , N , 7^ ■ (36) 
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Figure 2: The plot of EoS parameter w versus t. Here Hq — 0.1, m — I 



From equation ([35]), it is observed that the equation of state parameter w is time dependent, it can be function 
of redshift z or scale factor a as well (as already discussed in previous Section 1). 



So, if the present work is compared with experimental results (Knop et al. 2003; Tegmark et al. 2004b 
Hinshaw et al. 2009; Komatsu et al. 2009), then one can conclude that the limit of cj provided by equation ([35] 
may accommodated with the acceptable range of EoS parameter. Also it is observed that a.t t — tc, u! vanishes 
where tc is a critical time given by the following relation 



27 3n 



4(m + 2)^ \ tcj {m + 2)t^ 

Thus, for this particular time, our model represents a dusty universe. We also note that the earlier real matter 
at t < tc, where a; > later on at t > tc, where uj < converted to the dark energy dominated phase of universe. 

From Eq. (|34l) . we note that energy density of the fluid p{t) is a decreasing function of time and p > when 

Figure 1 is the plot of energy density of the fluid (p) versus time in accelerating mode of the universe. Here we 
observe that p is a positive decreasing function of time and it approaches to zero as t — cx). 

Figure 2 depicts the variation of EoS parameter (w) versus cosmic time (t) in evolution of the universe, as 
a representative case with appropriate choice of constants of integration and other physical parameters using 
reasonably well known situations (parameters are given in Figure caption). For m = 1, we obtain isotropic 
model which is studied here as a representative case. From Figure 2, we observed that at the initial time there 
is quintessence (w > —1) region and at late time it approaches to the cosmological constant (w = —1) scenario. 
This is a situation in early universe where quintessence dominated universe (Caldwell 2002) may be playing an 
important role for EoS parameter. Since w approaches to —1 for sufficiently large time, so its value is consistent 
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Figure 3: The plot of cosmological constant A versus t. Here £o = 0.1, m — 1 



with the range of all the three observations (Knop et al. 2003; Tegmark et al. 2004b; Hinshaw et al. 2009; 
Komatsu et al. 2009). 

In absence of any curvature, matter energy density flm and dark energy Q\ are related by the equation 

n^ + n^ = 1, (39) 

where ilm = 3^ and J^a = Thus, equation (p9)) reduces to 

SFT + sl^-'- <"»> 

Using Eqs. and in (PU)) . the cosmological constant is obtained as 

, 3 /5m + l\ n\2 3™_ 

A = --(^-^j(l + -) +M^V) (™... (41) 

From Eq. (|4ip . we observe that A is a decreasing function of time and it is always positive when 

nx2 3™ 44 / m + 2 



In general relativity, the Bianchi identities for the Einstein's tensor Gij and the vanishing covariant divergence 
of the energy momentum tensor together with imply that the cosmological term A is constant. In theories 
with a variable A-term, one either introduces new terms (involving scalar fields, for instance) in to the left hand 
side of the Einstein's field equations to cancel the non-zero divergence of Ag^ (Bergmann 1968; Wagoner 1970) 
or interprets A as a matter source and moves it to the right hand side of the field equations (Zcldovich 1968), in 
which case energy momentum conservation is understood to mean T**-' = 0, where T*j ~ Tij — (A/87rG).gij. It is 
here that the first assumption that leads to the cosmological constant problem is made. It is that the vacuum 
has a non-zero energy density. If such a vacuum energy density exists, Lorcntz invariance requires that it has 
the form (l),iy) = —{p)gfj.u- This allows to define an effective cosmological constant and a total effective vacuum 
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Figure 4: The plot of deceleration parameter q versus t 



energy density Ae// = A + 8ttG{p) or pvac = {p)+A/8nG. Note at this point that only the effective cosmological 
constant, Ag//, is observable, not A, so the latter quantity may be referred to as a 'bare'. The two approaches 
are of course equivalent for a given theory (Vishwakarma 2000). For detail discussions, the readers are advised 
to see the references (Carroll et al. 1992; Abdussattar and Vishwakarma 1996; Peebles and Ratra 2003; Sahni 
and Starobinsky 2000; Padmanabhan 2003, 2008). 

Figure 3 is the plot of cosmological constant A versus time t. We observe that cosmological parameter is 
decreasing function of time and it approaches a small positive value at late time (i.e. at present epoch). Recent 
cosmological observations (Perlmutter et al. 1998, 1999; Riess et al. 1998, 2004; Tonry et al. 2003) suggest 
the existence of a positive cosmological constant A with the magnitude A{Gh/c^) « 10""'^^'^. These observations 
on magnitude and red-shift of type la supernova suggest that our universe may be an accelerating one with 
induced cosmological density through the cosmological A-term. Thus, the nature of A in our derived DE model 
is supported by recent observations. 

Figure 4 is the plot of deceleration parameter q versus time t. From Figure 4, it is observed that q decreases 
very rapidly and reaches values — 1 then after it remains constant — 1 (like de Setter universe) . From this figure 
we observe that the DE model, for < n < 1.5, evolves from the matter dominated era to quintessence era and 
ultimately approaches to cosmological constant era where as for n > 1.5, the universe evolves from quintessence 
to cosmological constant era. It is worth mentioned here that for n < 1.5, transition of the universe takes place 
from the early decelerating phase to the recent accelerating phase where as for n > 1.5, the expansion of the 
universe is always accelerating. 

From these analysis we conclude that it is the choice of scale factor which makes the model inflationary at 
the early stages of the universe and radiation/matter dominance phase before the D.E. era. From Eq. (|29|) . 
we observe that when t — >■ 0, the expansion scalar becomes infinity which indicates the inflationary scenario. 
Also from Fig. 4, we observe that before t ^ 1, q > and this indicates radiation/matter dominance era of the 
universe. However, after t Ri 1, q < which indicates the DE dominated era. The solution in our model does 
not blow up at any given epoch for the choice of the ansatz (|20|) . Hence our derived model is physically acceptable. 
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The Cosmic Microwave Background (CMB) is also considered to be a major experimental evidence which 
support the present models of the observed universe and from this CMB observations several scientists found 
the signature of anisotropy. Based on these studies and observations, one may not preclude the possibility that 
our universe is anisotropic. We have already discussed this scenario in Introduction. 



700 



600 




m = 0.5 — - — m = I m = 2 



Figure 5: The plot of energy density p versus t. Here £o — 0.1 



5 STABILITY OF CORRESPONDING SOLUTIONS 

A rigorous analysis on the stability of the corresponding solutions can be done by invoking a perturbative 
approach. Perturbations of the fields of a gravitational system against the background evolutionary solution 
should be checked to ensure the stability of the exact or approximated background solution (Chen and Kao 
2001). Now we will study the stability of the background solution with respect to perturbations of the metric. 
Perturbations will be considered for all three expansion factors via 

ai ^ ast + Sai ^ OBiil + Sbi) (43) 

We will focus on the variables Sbi instead of Soi from now on for convenience. Therefore, the perturbations 
of the volume scale factor Vb — lif^iOi, directional Hubble factors 9i — ^ and the mean Hubble factor 

6 = X]i=3 ^ — W '^^^ ^® shown to be 

i i i 

One can show that the metric perturbations Sbi^ to the hnear order in 5bi^ obey the fohowing equations 

Y,Sb^ + 2J2SB^5h = (45) 

i 

Vr ■ ■ 

Sb, + -^Sb, + J2^^l^B^ = (46) 
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Figure 6: The plot of EoS parameter iv versus t. Here £o = 0.1 



From above three equations, we can easily find 



6b, + = 0, 

VB 

where Vb is the background volume scale factor. In our case, Vb is given by 

using above equation in equation (6) and after integration we get 

4n„-|tT" ^ ^ 3 13 



Shi = Cit"3"e"3*WittakerM ~-n,--n+ -, -i , 



(47) 



(48) 



(49) 



(50) 



where Ci is an integration constant. Therefore, the "actual" fiuctuations for each expansion factor 5ai — ub 
is given by 

3 3 13 



Saj Cit 4 e 4 WittakerM — n, n H — ,-t 

'4 4 2 2 



(51) 



From above equation we see that for n >> 1, 5ai approaches zero. Consequently, the background solution is 
stable against the perturbation of the graviton field. 



6 OTHER DARK ENERGY MODEL 

Now we take the following ansatz for the scale factor, where the increase in terms of time evolution is 

1 



a{t) 



t 



(52) 



12 




m = 0.5 



m = I 



m - 



Figure 7: The plot of cosmological constant A versus t. Here £o = 0.1 



By the above choice of scale factor yields a time dependent deceleration parameter and the corresponding 
solutions are stable. The motivations for selecting such type of scale factors for finding solutions are already 
described in Sect. 3. We define the deceleration parameter q as usual, 



aa 



Using ([521) into we find 



q = -2 



V2t3 + 1^ 

Using equations (IT6l) . (|24)) and ((52)) in (fTTj) . we obtain the expressions for metric functions as follows 



B{t) = £4, 



1 „\ 
1 



where, £4 = I ("'+^' , £5 = ££4, £q = £'^ and I is an integrating constant. 
Hence the model ([T]) reduces to 

6m 

1 



(53) 
(54) 

(55) 

(56) 
(57) 



ds' = -de + £2 f - - + dx^+£l{-- + A dy^ 



t 



t 



+ tl[--^+tA dz^ 



(58) 
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Figure 8: The plot of deceleration parameter q versus t 



The expressions for the Hubble parameter (i?), scalar of expansion (0), shear scalar (cr), spatial volume {V) and 
the average anisotropy parameter {Am) for the model (|58|) are given by 



3 f2t^ + 1 
iH^- \ - 



= 3 



1\ (2i3 + l) 



V 



Am, — 2 



n 2 



3 



TO — 1 



(59) 



(60) 



(61) 



(62) 



From Eq. ([5^ . we observe that when t ^ 0, 6* — ^ oo and this indicates the inflationary scenario at early stages 

2 

of the universe. Since ^ 7^ for all values of to except for to = 1, hence the model is anisotropic except for 
TO = 1. The the dynamics of the mean anisotropy parameter depends on the value of to. The mean anisotropic 
parameter is constant. We observed that when to — —2, Am 00 and for to = 1, Am = 0. Thus, the observed 
isotropy of the universe can be achieved in phantom model (see, Figure 6). 



The energy density of the fluid can be find by using Eqs. ([55l 



9(2m+ 1) {2t^ + 1)2 
(to + 2)2 (t3_ 1)2^2 



(157|) in (IT9|) 

6 m 
(m + 2) 



(63) 



where h = ji- Using Eqs. ([55]) . ([57]) and in pT)) . the EoS parameter uj is obtained as 



27 {2t^ + lf 6 (2t''+8f^-l) 

(m+2)2 ((3-1)2(2 (m+2) (t3-l)2t2 



1 t2\-(^;r+2) , 9(2m+l) (2/Hl) 



(64) 
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Using Eqs. (1551) . ([57]) . (|63p and ([M)) in ([T5)) . the skewness parameters S (or 7) (i.e. deviations from uj along 
y-axis and z-axis) are computed as 



(65) 



So, if the present work is compared with experimental results (Knop et al. 2003; Tegmark et al. 2004b; Hinshaw 
et al. 2009; Komatsu et al. 2009), then one can conclude that the limit of ui provided by Eq. ([M)) may 
accommodated with the acceptable range of EoS parameter. Also it is observed that at t = tc, uj vanishes, where 
tc is a critical time given by the following relation 

27 (2.? + I? 6 + - 1) , ^f_l_^ ^ ^ ^^^^ 



(m + 2)2 (^3 _ 1)2^2 + 2) [tl - 1)2^2 ^' y 

Thus, for this particular time, our model represents a dusty universe. We also note that the earlier real matter 
at t < tc^ where a; > later on at t > tc, where uj < converted to the dark energy dominated phase of universe. 

From Eq. (j63p . we note that energy density of the fluid p{t) is a decreasing function of time and p > when 

(2t^ + l)^ f 1 , ,,\™^ 4(m + 2)2 

(^3- 1)2^2 \ J - 9(2m + l) • ^ ' 

Figure 5 is the plot of energy density of the fluid (p) versus time t. Here we observe that p is a positive decreasing 
function of time and it approaches to zero as i — > 00. 

Figure 6 depicts the variation of EoS parameter (w) versus cosmic time (t) in evolution of the universe, as 
a representative case with appropriate choice of constants of integration and other physical parameters using 
reasonably well known situations (parameters are given in Figure caption). From Figure 6, we observe as follows: 

(i) for m < 0.5, the evolution of the universe starts from quintessence era {uj > —1) and approaches to 
phantom region (uj < —1). 

(ii) for 1 < m < 2, the universe evolves from phantom region (a; < —1), then crosses PDL and ultimately 
approaches to quintessence region (uj > —1). 

(iii) for 2 < m < 3, the evolution of the universe commence from phantom region (uj < —1), then crosses 
PDL and then skip over to non-dark region. 

(iv) for 3 < m, the evolution of the universe begins from quintessence era {uj > —1) and ultimately pass over 
to non-dark region. 

(v) for TO = 1, we get uj ~ —0.65 which is consistent with SNe la data —1.67 < uj < —0.62 (Knop et al. 2003). 

(vi) for TO = 0.5, we get uj ~ —1.1 which is reproducible with current observational realm (Knop et al. 2003; 
Tegmark et al. 2004b; Hinshaw et al. 2009; Komatsu et al. 2009). 

Using Eqs. (|59p and in (PH)) . the cosmological constant is obtained as 



-6n 



3(4to2 + 10to+ 13) (2^3 + 1)2 ^ ^ f 1 ^ ,2\'^^ 



^ (m + 2)2 {t^ - l)2t2 +^"[ 

From Eq. (j68p . we observe that A is a decreasing function of time and it is always positive when 

(2.3 + 1)2 ^ ™ ^ _^Jo(- + 2)- ^ , (69) 



(.3 - l)2t2 \ t J 3(4to2 + 10m + 13) ' 

Figure 7 is the plot of cosmological constant A versus time t. It is observed that in all cases cosmological param- 
eter is decreasing function of time and it approaches a small positive value at late time (i.e. at present epoch). 
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Thus, the nature of A in this derived DE model is also in good agreement with recent observations (Perlmutter 
et al. 1998, 1999; Riess et al. 1998, 2004; Tonry et al. 2003). 



Figure 8 is the plot of deceleration parameter q versus time t. From the figure we observe that the expansion 
of the universe starts from accelerating phase and the rate of expansion decreases with time and it stops and 
again starts accelerating to approach —0.5 which is very close to the value (~ —0.7) predicted by the observations 
(Riess et al. 2004; Virey et al. 2005). 

A convenient method to describe models close to A CDM is based on the cosmic jerk parameter j, a dimcn- 
sionless third derivative of the scale factor with respect to the cosmic time (Chiba and Nakamura 1998; Sahni 
2002; Blandford et al. 2004; Visser 2004, 2005). A deceleration-to-acceleration transition occurs for models with 
a positive value of jo and negative qo. Flat A CDM models have a constant jerk j — I- The jerk parameter in 
this case is obtained as In this case, we obtain the jerk parameter as 

2t'"^ + 2*4 - - t - 2 

This value is consistent with observational value j ~ 2.16 obtained from the combination of three kinematical 
data sets: the gold sample of type la supernovae (Riess et al. 2004), the SNIa data from the SNLS project 
(Astier et al. 2006), and the X-ray galaxy cluster distance measurements (Rapetti et al. 2007) for t = 1.50. 



7 STABILITY OF CORRESPONDING SOLUTIONS 

The method to study the stability of the background solution with respect to perturbations of the metric is 
already given in Sect. 5. From Eqs. (PSI) — (|T7)) . we can easily derive 

5h + = 0, (71) 

VB 

where Vb is the background volume scale factor. In our case, Vb is given by 

Vs^t^ (72) 

Using above equation in Eq. (j7ip and after integration we get 

5h = c,t-\ (73) 

where Ci is an integration constant. Therefore, the "actual" fluctuations for each expansion factor 5ai = asiSbi 
is given by 

5a, -> c,t-\ (74) 

where aBi — > ■ From above equation it is obvious that 5ai approaches zero as i — >■ oo. Consequently, the 
background solution is stable against the perturbation of the graviton field. 



8 CONCLUDING REMARKS 

A new class of anisotropic B-T^/o DE models with variable EoS parameter lo has been investigated by using time 
dependent deceleration parameter. In literature it is plebeian to practice a constant deceleration parameter. Now 
for a Universe which was deceleration in past and accelerating at present epoch, the DP must show signature 
flipping as already discussed in Section 2. Therefore our consideration of DP to be variable is physically justified. 

The main features of the models are as follows: 



• DE models present the dynamics of EoS parameter w provided by Eqs. (1351) and (|55p whose range are in 
good agreement with the acceptable range by the recent observations (Knop et al. 2003; Tegmark et al. 
2004b; Hinshaw et al. 2009; Komatsu et al. 2009). 
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• It can be easily seen that in both DE models, the mean anisotropic parameter vanishes at m = 1. Thus, 
our both anisotropic models approach to isotropy at m = 1. 

• We obtain cosmological constant dominated universe, quintessence and phantom fluid dominated universe 
(Chevallier and Polarski 2001), representing the different phases of the universe through-out the evolving 
process for different cosmic times. These fits suggest that u) > —1 for a long (quintessence-like) period in 
the past, and at the same time they suggest that the universe has just entered a phantom phase oj < —1 
near our present. 

• Unlike Robertson- Walker (RW) metric Bianchi type metrics can admit a DE that wields an anisotropic EoS 
parameter according to the characteristics. Therefore, one can not rule out the possibility of anisotropic 
nature of DE in the frame-work of B-VIq space-time. 

• In first case, the observed isotropy of the universe can be achieved in cosmological constant model (see. 
Figure 2) whereas in second case, the observed isotropy of the universe can be achieved in phantom model 
(see. Figure 6). Thus, Bianchi type-VIo models which remain anisotropic are of preferably academical 
interest. 

• Our DE models are of great importance in the sense that the nature of decaying vacuum energy density 
A{t) is supported by recent cosmological observations (Perlmutter et al. 1998, 1999; Riess et al. 1998, 
2004; Tonry et al. 2003). 

• Though there are many suspects (candidates) such as cosmological constant, vacuum energy, scalar field, 
brane world, cosmological nuclear-energy, etc. as reported in the vast literature for DE, the proposed 
models in this paper favour EoS parameter as a possible suspect for the DE. 

• The cosmic jerk parameter in our derived models is also found to be in good agreement with the recent 
data of astrophysical observations namely the gold sample of type la supernovae (Riess et al. 2004), the 
SNIa data from the SNLS project (Astier et al. 2006), and the X-ray galaxy cluster distance measurements 
(Rapetti et al. 2007). 

• For different choice of n, wc can generate a class of DE models in Bianchi typc-VIo space-time. It is observed 
that such DE models are also in good harmony with current observations. Our study is continued and we 
shall generate some other interesting physically viable models for other values on n. 

• Our corresponding solutions have inflationary scenario at the early stages of the universe and also radia- 
tion/matter era before DE era. 

• Our corresponding solutions are physically acceptable and the solutions are stable. 

Thus, the solutions demonstrated in this paper may be useful for better understanding of the characteristic of 
anisotropic DE in the evolution of the universe within the framework of Bianchi type-VIo space-time. 
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